Characterization and Corrections of Relative Humidity Measurement from Meteomodem M10 Radiosondes at Midlatitude Stations

Measurement of water vapor or humidity in the atmosphere is fundamental for many applications. Relative humidity measurements with a capacitive sensor in radiosondes are affected by several factors that need to be assessed and corrected. This work aims to address corrections for the main effects for the MeteomodemM10 radiosonde as a step to meet the Global Climate Observing System (GCOS) Reference Upper-Air Network (GRUAN) requirements. The considered corrections are 1) the calibration correction; 2) a slow regime due to the slow diffusion of molecules through the sensor, especially at very high and very low relative humidity conditions; 3) the relative humidity sensor dependence on the gradient of temperature; and 4) the time lag at cold temperatures, which affects measurements in regions of strong relative humidity gradients. These corrections were tested for 26 nighttime and 25 daytime radiosondes in two midlatitude locations for which both MeteomodemM10 and Vaisala RS92 measurements were available. The results show that, after correcting for the four effects, M10 relative humidity measurements are, on average, consistent with the Vaisala RS92 relative humidity values within 2% RH at all altitudes for the nighttime launches (against 6% RH before the correction) and within 5% RH at all altitudes for the daytime launches (against 9% RH before the correction).


Introduction
Measurement of water vapor or humidity in the atmosphere is fundamental for many applications (Bojinski et al. 2014). Humidity measurement at higher altitude is still among the most difficult problems in basic meteorology (Elliott et al. 2002): the measurement depends on the exchange of water molecules between the sensor and the air, and during the sounding water vapor concentrations decrease by several orders of magnitude between the surface and the tropopause. Accurate humidity measurements from the ground to the top of the troposphere and lower stratosphere are necessary to prevent spurious drying or moistening of the atmosphere simulated in numerical weather prediction models (Seidel et al. 2009). Such measurements are also necessary to better understand the role of water vapor in climate feedback, to study the life cycle of anthropogenic clouds such as contrails, and to quantify the state of deliquescence of atmospheric aerosols, just to name a few among many other fields of investigation. Moreover, vertical profiles of temperature and humidity can be used to calibrate (error and uncertainty) different types of satellites with active and passive remote sensing instruments (Zou et al. 2013).
Today, the most widely used atmospheric humidity measurements rely on capacitive sensors that are based on a dielectric polymer deposited between two electrodes. The change in the dielectric constant of the polymer is directly proportional to the change in relative humidity, especially in the 5%-95% relative humidity range (Balagurov et al. 1998). Some instruments also use sensors that are based on the impedance of conductive polymers, which shows an exponential relationship to relative humidity. Capacitive sensors are typically preferred because of a lower sensitivity to temperature of the absolute calibration (Matsuguchi et al. 1998).
The relative humidity measurement accuracy is generally considered to be 2%-5% RH, and this variability is due to 1) errors in absolute calibration; 2) efficiency of the diffusion of water vapor molecules inside the sensor; 3) capacitance sensor response time, which ranges from seconds to minutes depending on temperature; and 4) bias due to spurious heating or cooling of the capacitance sensor relative to ambient air.
All these effects on relative humidity measurements must be well understood and accounted for to characterize measurement accuracy. Each error and its correction has an associated uncertainty that have been carried out and published for several Vaisala, Inc., radiosonde sensors (RS80: Miloshevich et al. 2001Miloshevich et al. , 2004RS92: Vömel et al. 2007;Dirksen et al. 2014), and for the Meisei Electric Co., Ltd., radiosondes (RS-11G and IMS-100: Kizu et al. 2018).
In the framework of the Global Climate Observing System (GCOS) Reference Upper-Air Network (GRUAN), Dirksen et al. (2014) propose methods to correct Vaisala RS92 relative humidity measurements and to provide measurement uncertainty estimates. Such methods are applied to improve data quality at all GRUAN stations that use RS92 radiosondes. Similar investigations should be carried out on other relative humidity sensors of radiosondes used in the GRUAN network. At least two stations of the GRUAN network (Paris, France and Réunion in the southern Indian Ocean) use Meteomodem Co. M10 radiosondes for regular temperature and relative humidity profiling.
The aim of this paper is to identify error sources that can lead to dry or wet bias and uncertainties. This paper presents a set of corrections that aim at reducing this bias and uncertainties using intercomparison measurement with the RS92. Section 2 presents the dataset used in the study and provides a description of each of the correction methods. Section 3 discusses the impact of each correction on M10 measurements and describes the general result when compared with RS92 measurements. Section 4 presents the sensitivity of measurement bias and standard deviations to different values of correction parameters.

Observational dataset
To test correction methods and evaluate their effects on measurement accuracies, we looked for a dataset that contains a significant number of both nighttime and daytime dual-radiosonde profiles that include M10 radiosondes and another well characterized radiosonde that can serve as reference, both under the same balloon. This well characterized radiosonde is here defined as a radiosonde for which there exist published studies on its accuracy and on its measurement uncertainties. This is the case for Vaisala RS92 radiosondes, which have been extensively tested (Dirksen et al. 2014). For this study, the data from two field experiments, which are presented in the following subsections, are considered.

a. Meteomodem M10 radiosonde
The M10 is one of the radiosonde commercialized by Meteomodem. The M10 was released in 2010, and it is fully compatible with the SR10 receiver system and the Meteomodem software (Fig. 1). Its dimensions are 95 mm 3 95 mm 3 88 mm for 150 g with batteries. The M10 is composed of a capacitive humidity sensor covered by an innovative metal coated shield that allows good ventilation while protecting it from direct radiation and from water droplet freezing on the sensor, a temperature sensor that measures the air temperature and positioned at the very end of the sensor boom, and a GPS sensor whose coordinates are measured and from which the pressure, the vertical velocity, the wind speed and direction are derived. The capacitive humidity sensor is composed of three primary components: a basic layer that acts as an electrode; a dielectric material, whose characteristics are a function of relative humidity; and a fast response porous electrode that acts as the

858
second electrode of the capacitor. A second thermistor is located under the protective shield close to the humidity sensor in order to have an approximative measurement of the temperature of the capacitive humidity sensor. The UPSI France Company is the subcontractor for this capacitive humidity sensor and these sensors are made specifically and exclusively for Meteomodem. Figure 2 shows the localization of all active sites using Meteomodem radiosonde technology in 28 countries, colors are relative to the year of the station deployment.  Bock et al. (2013) found a dry bias for M10 observations of 26.5% RH in the altitude range of 1-12 km and increasing to 210% RH at 10 km when compared with Meteolabor Snow White.
For that study, only calibration-related corrections were applied to the M10 measurements, as the M10 radiosondes had been released shortly prior to this campaign. The experience acquired with M10 afterward during other campaigns and the willingness to use M10 radiosondes at several measurement stations of the GRUAN network, led to the development of the corrections that are presented in the section 3.
Seven Meteolabor Snow White (night models; 026 and 059 versions) soundings and 20 RS92 soundings were performed during the campaign (Table 1). In this study, we use these 20 soundings to apply the corrections to the M10 RH and perform comparisons with the RS92 RH profiles. For these soundings, the Vaisala DigiCORA III software, version 3.64, was used. No additional correction to the default corrections were applied to the RS92 data (Bock et al. 2013). According to the Vaisala data continuity website (https://www.vaisala.com/en/press-releases/ 2009-12/vaisala-launches-website-ensure-data-continuityradiosounding-observations), the time lag, the daytime correction for solar heating, and the nighttime correction for radiative cooling started to be applied in December 2010 for standard RH data products recorded with DigiCORA 3.64 software. Because all measurements were performed at night during the OHP campaign, the only correction applied to RS92 humidity was the time-lag correction, the solar radiation algorithm having no impact on nighttime humidity (Jauhiainen et al. 2011). The accuracy of all of these corrections applied to the RS92 data used in this study in comparison with the GRUAN-corrected RS92 data is discussed in Yu et al. (2015).

c. SIRTA field experiment
The Instrumented Site for Atmospheric Remote Sensing Research (SIRTA) is a French national observatory dedicated to studying clouds, aerosols, dynamics and thermodynamics in the boundary layer and the free troposphere. The SIRTA observatory is a midlatitude site (48.78N, 2.28E) located in a semiurban area, on the Saclay plateau 25 km south of Paris (Haeffelin et al. 2005).
Three intensive observational periods (IOPs) dedicated to comparisons of vertical profiles of humidity were carried out between October 2016 and August 2018 with 29 M10-RS92 dual flights (see more details in Table 1). In this study, both IR2010 and Eoscan software were used to record the data measured and sent by the M10 radiosondes. IR2010 is the former software, which is still in used at several sites, whereas Eoscan is the new software from Meteomodem. For this study, we used the Eoscan and IR2010 programs solely to produce the raw data before any posttreatment. DigiCORA III software was used for the monitoring of the RS92 radiosondes and we allowed the standard corrections to be applied for the data considered in our study. The different corrections presented in section 3 are applied to the raw M10 data.
Standard ground checks were applied for M10 (ambient temperature and humidity chamber) and RS92 (temperature and 0% RH humidity chamber). For all of the figures, the altitude corresponds to the altitude above mean sea level. Figure 3 shows a schematic overview of the correction methodology to convert the raw value of relative humidity to an optimized value taking into account the four corrections presented in the paper. Input parameters and related equations are indicated.

Correction methods
Meteomodem has developed one transfer function and performs an absolute calibration of the capacitance sensor that allows the conversion of the measured frequency into a relative humidity. This calibration can be verified a posteriori by measuring the relative humidity in conditions where the relative humidity is known. The sensor response time to a step change is highly dependent on temperature following a power law, yielding time constants that exceed 1 min when the temperature is below 2558C. Miloshevich et al. (2004) propose a response-time correction that we test and implement. It is observed that the capacitance relative humidity sensors are also affected by an additional slow-regime response time that has not been discussed in the literature. We present the concept and propose a correction method. Last, with regard to relative humidity bias due to spurious temperature of the sensor, often referred to as ''radiation dry bias,'' several authors (e.g., Vömel et al. 2007;Miloshevich et al. 2009;Dirksen et al. 2014) propose correction methods that depend on temperature and relative humidity or pressure. These methods are used and tested in the analysis of the data (section 3c).

a. Absolute calibration correction
The general process of calibration for the M10 radiosondes concerning the relative humidity sensor (RHS) follows two separate procedures. In the factory, each radiosonde . The SCS is equipped to provide up to five predefined humidities in the 0.5%-85% RH range at 238C with a accuracy better than 0.6% RH. It is also equipped with an air conditioned chamber allowing to test up to 20 different measurement ranges with a 2s tolerance smaller that 2.1% RH (Accreditation Directory SCS0065). The M10 relative humidity sensor is calibrated for this 0.5%-85% RH range. According to Rotronics, the HC2S is traceable to an international standard through a documented and unbroken chain of calibrations. The HC2A-S used in the factory calibration process is compared four times per year with one other HC2A-S devices to verify the absence of measurement deviation. An additional verification procedure is carried out during the Meteomodem ground check before each launch. The ground check humidity sensor is calibrated with respect to a reference standard directly in laboratory conditions. For this study, concerning the radiosondes launched during the OHP field experiment in 2011, a supplementary correction for calibration has been performed on each RH sounding record, which consists in making the assumption that the RH values converge to 0% at 30-km altitude, as expected. Each radiosonde has been verified individually to be sure that there is no contamination of the RH sensor by frozen water droplets. This correction had significant values for the case of the M10 radiosondes used in OHP campaign, which happened to be the prototype series used by Meteomodem for this type of radiosondes (M10 replaced the M2K2 radiosondes after that campaign) and for which issues during the early phases of the M10 radiosonde industrial production were encountered. This supplementary correction was not needed for the radiosondes used for the launches at SIRTA. FIG. 3. Schematic overview of the correction method: TS is the temperature of the humidity sensor, U 0 is the raw humidity value, U 4 is the final humidity value after all of the corrections, T 0 is the raw temperature value, and T 1 is the air temperature after the absolute calibration if needed.
MAY 2020 D U P O N T E T A L .

b. Slow-regime correction
A new method has been developed to take into account the effect inspired by the layered model of a capacitive humidity sensor proposed by Wildmann et al. (2014). Our method separates the response of the sensor to a step change in RH into a combination of two layers: a fast-regime layer and a slow-regime layer. The fast-regime layer corresponds to a fast change in humidity by more than 90% of the step change within a few seconds. The slow-regime layer represents a slow evolution of the measured humidity in the range of minutes for about 10% of the step change. To correct the measured humidity, a fraction of the difference between the measured humidity and the slow-regime layer is added to the measured humidity at each time step. We can define the corrected relative humidity U 2 as follows: where U slow is the relative humidity in the slow-regime layer, U 1 is the measured and uncorrected relative humidity, U 2 is the corrected relative humidity, m is the fraction of the slow regime, and i is the index for a given position; where t (in seconds) is the time constant of the slowregime part of the sensor. We have tested many combinations of t and m, and we show some examples in Fig. 4 [(120 and 0.2); (500 and 0.1); (1000 and 0.05)]. The left panel of Fig. 4 shows U slow for different values of t (120, 500, and 1000 s) with a constant m 5 0.1. The right panel of Fig. 4 shows U slow for different values of m (0.2, 0.1, and 0.05) with a constant t 5 500 s, for a particular example profile.
As expected, the higher the variation of relative humidity is, the more important the correction will be, especially for a small value of t, that is, relatively fast response of the slow-regime layer, emphasized by an important m. Slow-regime correction will be effective just after crossing a cloud under a dry air mass and, to a relatively smaller extent, after crossing the tropopause (see Fig. 4 at 1.8 and 12 km, respectively).

c. Temperature-gradient correction
For a good measurement, the relative humidity should be measured when the capacitive sensor temperature and the air temperature are the same. However, multiple factors may hinder or contribute to break this equilibrium. During daytime measurements, as for the ones at the SIRTA observatory, the solar radiation can heat the sensor and create a dry bias in the humidity measurements, whereas during nighttime measurements, as at the OHP, the sensor is exposed to infrared radiation from Earth and from the surrounding water vapor and clouds. These radiative effects are modulated both in daytime and nighttime periods with the heat transferred by conduction to the sensor, from the radiosonde to the sensor boom and with the ventilation leading to a cooling effect.
Each of these factors has a variable influence on the sensor heat budget and, therefore, on its temperature. Indeed, radiative effects are very variable along the vertical due to the role of clouds on direct solar beam or on infrared cooling in relation with water vapor and hydrometeor spatial distribution. The ventilation of the sensor depends on the ascent speed, the distance between the radiosonde and the balloon and the induced pendulum motion, the presence of other radiosondes and their configuration under the balloon.
To evaluate the heating of the RHS, the approach undertaken for the M10 radiosondes was to take advantage of an available temperature measurement close to the RHS (hereinafter T S ) and to compare it with the measured temperature of the air T 1 . Temperature T S is measured on the sensor boom, 2 mm away from the RHS edge. Since the temperature sensor is not directly located on or inside the RHS, we can directly consider 1) the instantaneous value of the temperature difference between the air and the RHS (T 1 and T s measurements) and 2) a fitted model based on T 1 and T S . Figure 5 shows an example of this temperature difference profile (T S 2 T 1 , denoted as DT) at the SIRTA (Fig. 5, left panel) and OHP (Fig. 5, right panel) sites. The measured DT can be fitted with a polynomial fit as follows: where c i coefficients are optimized by the least squares method, and z is the altitude above ground level in kilometers.
The c i coefficients can be adjusted for each flight and for all the daytime or nighttime periods. The main advantage for this daily modeling is to have a better accuracy (see Fig. 5, blue lines), but the drawback is that we are obliged to await the end of the flight to process the DT. The main advantage of the climatological modeling is that we are ready to process DT at the beginning of the launch but the main drawback is that the uncertainty is bigger. In Fig. 5, the black line corresponds to the measured value; red, green and blue lines correspond to the different models (daily or climatological).  Table 2 shows 1) the optimal parametrical coefficients for the daytime and nighttime climatological models obtained for the Trappes, France, radiosoundings between 2012 and 2016 (a total of 1800 radiosoundings at 1200 and 0000 UT), and 2) the coefficients for the two selected dates at SIRTA and OHP. For the climatological models, we have gridded all the profiles between the surface and 12 km of altitude with a vertical resolution of 50 m, and we optimized the parametrical coefficients by the least squares method. To account for the instantaneous effect of clouds, solar fluxes and convection effect that are shown on Fig. 5 (instantaneous variability reaching more than 108C and a relative humidity correction exceeding more than 10%), we consider in this study the instantaneous value of the temperature difference between the air and the RHS rather than daily or annual parametrical model. For the SIRTA radiosounding during daytime, the difference of temperature between the RH capacitive sensor and the air temperature sensor ranges from 08 to more than 68C, at the surface and at 12 km, respectively. The solar radiation induces a difference that is strictly positive, with a humidity sensor warmer than the air, leading to a dry bias that needs to be corrected. To minimize the effect, the shield cover shades the humidity sensor. Because of its design, the humidity sensor is black, so it is expected that the influence of infrared radiation will be high, unlike air temperature sensor, which is covered by a metallic coating. The protection therefore shades the RH sensor to limit the effect of radiation. Its shape and design have been produced by Meteomodem.
For the OHP radiosounding during nighttime, DT ranges from 218 to 128C, at the surface and around 12 km, respectively. The infrared radiative cooling associated with the dynamics around the M10 radiosonde induces alternatively a warming or a cooling of the humidity sensor. We note that the daily fit has a very good agreement with the measured DT whereas the climatological model can have discrepancies reaching 18 and 0.58C for daytime and nighttime periods, respectively.
The estimation of the relative humidity sensor heating allows us to correct the observed humidity using the measured difference between the temperature of the capacitive sensor and the true temperature of the air. The relative humidity is defined as the ratio of the water vapor partial pressure e and the water vapor saturated pressure e sat [Eq. (4), below]. The water vapor saturated pressure is a function of the temperature only; U 2(i) 5 e 2(i) e sat2(i) and U 3(i) 5 e 3(i) e sat3(i) , where the subscripts 2 and 3 correspond to two different locations on the sensor boom, and i is the index for a given position. The first corresponds to the location of the air temperature sensor (subscript 2), and the second corresponds to the location of the temperature sensor near the humidity sensor (subscript 3). Assuming that the water vapor partial pressure is uniform between these two locations (i.e., e 2 5 e 3 ), Eq. (4) leads to e 2(i) 5 e 3(i) $ U 3(i) 5 U 2(i) e sat2(i) e sat3(i) .
The Hyland and Wexler law (Hyland and Wexler 1983) leads to where T 2(i) is the air temperature sensor in kelvins, T S(i) is the temperature sensor near the humidity sensor in kelvins, A is equal to 20.580 3 10 4 , B is equal to 20.0486, C is equal to 10.418 3 10 24 , D is equal to 20.145 3 10 27 , E is equal to 16.546, and i is the index for a given position. We tested 10 different formulations of water vapor partial pressure over liquid water and decided to select the abovementioned Hyland and Wexler law (Hyland and Wexler 1983). For a 28C heating, which is a typical value of DT at 6 km for daytime periods, K is around 1.25. Likewise for a 48C heating, which is a typical value of DT at 10 km for daytime periods, K is around 1.5. The impact on the vertical profile of humidity can also be very significant, particularly for high altitude.

d. Sensor-response-time correction
The time constant of capacitance sensors to a step input change is typically of less than 1 s at 208C. The radiosonde ascent speed is typically 5 m s 21 ; hence, a 1-s time response would be sufficient to monitor the abrupt changes in relative humidity observed in the atmosphere at the top of the boundary layer or at cloud boundaries. However, the time constant increases significantly at colder temperatures, exceeding 30 s at temperatures below 2408C according to Meteomodem. Because the response time increases with decreasing temperature, sharp relative humidity changes in the upper troposphere and lower stratosphere are likely to be reported by the radiosonde at a later time than they truly occur. This effect must be corrected. Miloshevich et al. (2004) documented an approach to correct the response time for the Vaisala RS80 humidity sensor in which the response is an exponential function of time [see Eq. (3) in Miloshevich et al. 2004]. Dirksen et al. (2014) more recently developed a correction for the RS92 sensor in which the sensor response time of the humidity sensor is an exponential function of the temperature. Although the response-time effect is common to all capacitive sensors, the relative humidity sensor time-constant temperature dependences of two different sensors are likely not to be the same.
Considering the relative humidity sensor as a firstorder low-pass filter, the measured relative humidity can be expressed as where U 3 is the measured relative humidity, U 4 is the final relative humidity after all of the corrections, k is the gain factor, and i is the index for a given position. This equation can be mathematically rewritten, expressing U 4 as a function of the measured relative humidity slope: where t (T) is the temperature-dependent time constant of the sensor [see Eq. (10), below], N is the number of measurements used to derived the slope, and i is the index for a given altitude. This number N depends on the temperature and is defined as where a 1 5 20, g 1 5 4, and b 1 5 60; N ranges from N (T.08C) 5 5 s to N (T,2608C) 5 120 s to prevent extreme values at low altitude and in the upper troposphere.
Time constant values are provided by Meteomodem (from laboratory tests) for temperatures ranging 1208 to 2408C (see black circles in Fig. 6, left panel). They are not provided at lower temperatures, for which an empirical approach was used to derive the t (T) relation. For this, some radiosoundings performed by MeteoFrance at Trappes were used. The method retrieves t around the tropopause, where temperatures can be as low as 2708C at midlatitudes, by comparing the measurements during radiosonde ascent and descent at this altitude. The soundings were selected following three criteria: 1) nighttime soundings only, to avoid an additional effect from solar radiation, 2) a small stratospheric drift, to maximize the chances to probe the same temperature and humidity when passing upward and downward at the tropopause, and 3) a neat humidity drop at the tropopause level; for this, cases of deep convection were discarded to have stable conditions between upward and downward situations.
During both the ascent and descent of the radiosonde, the measured RH at the tropopause level was corrected following Eq. (10). The left panel of Fig. 6 shows values provided by Meteomodem and the values derived from the experimental method. The value of t [Eqs. (8) and (10)] is changed by iteration to minimize ''d,'' the Euclidean distance between the ascent and descent corrected curves around the tropopause.
To take into account these points and the ranges, a curve was derived using the following equation: where the values for the coefficients are a 5 90, g 5 4, and b 5 60. This method turns out to be a useful evaluation technique as an alternative to a retrieval of the RHS MAY 2020 D U P O N T E T A L .
response time in the laboratory. It is important to stress that the temperature dependence expression (Fig. 6, left panel) derived here is an evaluation of the RHS time-lag behavior with measured air temperature. Indeed, Meteomodem points at 2408 and 2208C are not perfectly superimposed to the curve. Neither are some of the experimental points. A compromise was then adopted to optimize the three parametrical coefficients to the complete dataset. Figure 6 (right panel) shows an example of the application of the time-lag correction on the RH measurements from 8 to 14 km of altitude for a M10 radiosounding. The black and red curves are, respectively, the RH before and after the time-lag correction. The correction is seen to enhance the slopes and the amplitudes of the peaks.

a. Sensitivity tests
We have presented in section 3 the method to optimize the M10 radiosonde profiling with three main corrections. The first correction, hereinafter ''slow-regime correction,'' depends on two parameters t and m. The second correction, hereinafter ''temperature-gradient correction,'' is related to the period during which the temperature offset is fitted (daily, climatological) or measured (instantaneous). The third correction, hereinafter ''sensor-response-time correction,'' is directly fitted with a previous dataset.
We estimate the impact of each successive corrections of M10 relative humidity (U 1 , U 2 , U 3 , and U 4 ) in comparison with RS92 relative humidity. In Tables 3 and 4, we report a global correction effect difference constructed as the absolute mean over all the flights and altitudes on an interpolated 50 m grid. The RS92 dataset has been used to develop the three main corrections on M10 radiosonde, and this dataset validates the new M10 dataset. Table 3 corresponds to the data collected at the SIRTA observatory with 25 radiosondes during the daytime periods. The initial difference in relative humidity between M10 and RS92 is around 5.4% RH with a standard deviation of 1.8% RH. The impact of t and m appears directly on U 2 , but the final impact on U 4 is very significant (U 4 ranges from 3.9% to 5.1% RH on average). The standard deviation is also divided by 2 between U 1 and U 4 . The best configuration corresponds to t 5 500 s and m 5 0.1 for a measured temperature gradient with a final mean difference of 3.9% 6 0.8%

866
RH. The daytime climatological model dedicated to fit with the temperature correction leads to a decrease in the standard deviation to 1.1% RH, and the average bias decreases from 5.4% to 3.9% RH. Table 4 is the same as Table 3 but for the OHP field experiment with 20 nighttime radiosoundings. The initial bias is equal to 3.2% RH with a standard deviation of 1.1% RH. The same configuration of parameters (t 5 500 s and m 5 0.1) is used to reach the minimum bias (1.7% RH) and standard deviation (0.6% RH), which both decrease by a factor 2. Figure 7 shows the vertical profile of slow-regime correction (left panels), temperature correction (center panels), and time-lag correction (right panels) for the 25 radiosoundings launched at the SIRTA observatory (Fig. 7, top panels) and for the 20 radiosoundings launched at the OHP (Fig. 7, bottom panels). The average value is plotted as a solid line, and the standard deviation is shown with dashed lines. The slow-regime correction depends on the difference between the humidity in the slow reservoir and the measured humidity. Its impact is more perceptible when RH changes are fast and important. The slow-regime correction subtracts a maximum of around 2.5% of RH at 5 km for the SIRTA profiles, and adds 2.5% of RH at 1 km and subtracts 1% of RH at 3 km for the OHP profiles. At higher altitudes, the slow-regime correction is smaller than 1% of RH. The curves in the center panels of Fig. 7 show the mean impact of the temperature-gradient correction. For the SIRTA field experiment (i.e., for daytime flights), between the ground and 5 km, this correction is around 13%-4%, meaning that on average the relative humidity sensor was warmer than the air at these altitudes with a relatively constant offset. Above 5 km, the mean impact of the correction increases and peaks at 10 km with 112% RH as a result of the solar heat flux associated with a high value of relative humidity (see section 4b), then decreases to 0% RH at 13 km, where the atmospheric relative humidity tends to 0% RH. For the OHP field experiment (i.e., for nighttime flights), this correction ranges from 0% RH at 5 km to 3% RH at 12 km or to 1% RH at 15 km in relation to the absence of solar heating.
The curves in the right panels of Fig. 7 show the impact of the time-lag correction. Up to 7 km, the correction is zero because the temperature is higher than 2408C, so that the sensor is still responding fast (see Fig. 6). Above this altitude, the correction can lead to an increase or decrease in RH, depending on each single profile slope (Fig. 6). The correction maximizes around 10-12 km for SIRTA and 11-13 km for OHP, corresponding to the fast relative humidity drop near the tropopause, with a maximum mean effect of 23% RH for SIRTA and 24% RH for OHP.
The blue line corresponds to the vertical profile of relative humidity difference between the M10 sensor before all of the corrections and the RS92 sensor. The cyan line corresponds to the vertical profile of relative humidity difference between the M10 sensor after all of the corrections and the RS92 sensor. For all of these corrections, the mean impact is processed after an interpolation to a 50-m vertical grid between the surface and the maximum altitude for which both radiosondes (M10 and RS92) have reported measurements. For the SIRTA field experiment, the initial (before the corrections) maximal difference is around 28% near the ground level, 13% RH at 5 km, and 26% RH at 10 km. After all of the corrections, this difference is around 24% near the ground level and reaches 15% at 7 km. For the OHP field experiment, the initial maximal difference is 23% RH at 1.5 km and reaches 26.5% RH at 11.5 km. After all of the corrections, the difference ranges between 22% and 12% RH. Figure 9 shows the mean relative humidity profiles and standard deviations for M10 (blue) and RS92 (red) interpolated on a 50-m vertical grid, for SIRTA ( Fig. 9, left panel) and OHP (Fig. 9, right panel) field experiments. The mean difference is smaller (with respect to each single flight) with a high stability of the difference (a low standard deviation relative to the first series of M10 data; Bock et al. 2013). The correction is much more efficient for nighttime profiles (OHP campaign) in relation with the temperature gradient that is very low for nighttime periods. For the daytime profile obtained during the SIRTA field experiment, the most important difference between M10 and RS92 (RH M10 . RH RS92 ) appears between 5 and 10 km, where a lot of cirrus clouds occur during the different field experiments. In some events, the temperature-gradient correction can have a too strong impact on the relative humidity profiles compared to the RS92 radiosonde.

Conclusions
A set of corrections for M10 RH has been developed and tested for two field experiments, one with daytime launches at the SIRTA observatory and one with nighttime launches at the OHP. The corrections are based on observed physical phenomena and postprocessing verification. The corrections applied on the dataset are the following: 1) a correction for calibration, applied to the raw data of the OHP campaign (necessary because of a calibration bias in the early series of M10 radiosondes), 2) a slow-regime correction, based on the observed behavior of the sensor when responding to a step change in RH, FIG. 7. Vertical profile of (left) the slow-regime correction, (center) the temperature correction, and (right) the time-lag correction for (top) the 25 radiosoundings launched at the SIRTA observatory for the daytime period and (bottom) the 20 radiosoundings launched at the OHP Observatory for the nighttime period. The average value is plotted as a thick line, and the standard deviation is shown with dashed lines.
3) a temperature-gradient correction, which is based on the estimation of the difference between the RH sensor temperature and the actual air temperature, and 4) a time-lag correction, which improves the RH sensor response at very low temperatures.
The corrections have been applied to the whole M10 dataset and the impact of each correction on the M10 relative humidity profile has been quantified. As mentioned by Bock et al. (2013), the M10 used during the OHP campaign was a very recent product at the time. These radiosondes were actually the very first series that had been produced. The calibration issues encountered then, because of the early production stage, have a priori no link with the physical nature of the M10 relative humidity sensor.
The three remaining corrections are based on physical phenomena. The time-lag correction and the slowregime correction have comparable amplitude for day and night periods, but the time-lag impact is much more significant at high altitude (above 10 km) whereas the slow regime affects the relative humidity in the lower to middle atmosphere (below 7 km).
The time lag affects all capacitive relative humidity sensors and is due to the temperature dependence of the sensor response time at temperatures lower than approximately 2408C. Here, an original method using ascent and descent profiles has been used to derive an estimation of the dependency of the sensor response time to temperature.
The temperature-gradient correction uses a temperature measurement that is done near (2 mm away) the RHS. The magnitude of the correction factor depends on the measured temperature difference between the air and the sensor. It allows a correction of the measured RH, which can be positive or negative depending on the sign of the difference.
The observed slow reaction of the RH sensor when submitted to a relative humidity step change has motivated the development of a slow-response correction. This correction is based on the assumption that the propagation of water vapor molecules inside the relative humidity sensor can be explained with the diffusion theory equations.
Meteomodem developed a transfer function and performs an absolute calibration of the capacitance sensor that allows the conversion of the measured frequency into a relative humidity. This calibration can be verified a posteriori by measuring the relative humidity in conditions where the relative humidity is known. The sensor response time to a step change is highly dependent on temperature following a power law, yielding time constants that exceed 1 min when FIG. 8. Average contribution of each correction (slow regime, black; temperature, green; time lag, red) to the vertical profile of relative humidity (raw data in blue, and final result in cyan) for (left) the 25 radiosoundings launched at SIRTA for the daytime period and (right) the 20 radiosoundings launched at OHP for the nighttime period. the temperature is below 2558C. Miloshevich et al. (2004) propose a response-time correction that we test and implement. It is observed that the capacitance relative humidity sensors are also affected by another effect that can be explain by an additional slow-regime response time that has not been discussed in the literature. We present the concept and propose a correction method. With regard to relative humidity bias due to spurious temperature of the sensor, often referred to ''radiation dry bias,'' several authors (e.g., Vömel et al. 2007;Miloshevich et al. 2009;Dirksen et al. 2014) propose correction methods that depend on temperature and relative humidity or pressure. We test and implement such methods. The average RH profiles measured by M10 and RS92 have been compared after each correction has been applied to the M10 observations. The remaining difference between the two sensors is lower than 1.7% RH and 3.9% RH for the OHP and SIRTA field experiments, respectively.
In a future work, we will use these three main corrections to evaluate the Meteomodem M10 final products with 1) the GRUAN-corrected RS92 dataset, 2) a Snow White radiosonde, and 3) a cryogenic frostpoint hygrometer (CFH) radiosonde for more than 10 dual flights conducted at Réunion Island between 2017 and 2019.